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In an recent paper™ Regina, Rastelli and Tassi inves- 
tigated the critical properties of a two dimensional Ising 
model with exchange and dipolar interactions described 
by the Hamiltonian 



T~l = — J ^ ' <Ji<Jj 



<i,3> 



^3 



(1) 



for different values of the ratio J/g using Monte Carlo 
(MC) simulations and finite size scaling for different sys- 
tem sizes. Hamiltonian |T]) describes an ultrathin metal- 
on-metal magnetic film with perpendicular anisotropy in 
the monolayer limit 2 -. Almost simultaneously we pub- 
lished a paper— about the same subject with a very sim- 
ilar analysis for a different set of values of the ratio J/g. 
Our detailed numerical analysis showed evidences of an 
intermediate disordered phase, which was not detected 
by Rastelli analysis. The particular characteristics of 
that intermediate phase and its associated phase tran- 
sitions can introduce a strong bias in the results if not 
properly taken into account by the numerical procedure. 
In this comment we show that the numerical calculation 
protocol used by Rastelli et al is unable to detect the 
presence of the intermediate phase leading to spurious 
results. Therefore some of their conclusions must be re- 
vised. 

The main point concerns the results of Rastelli et al 
for J/g = 3.4 and J/g = 5 compared with our results 
in the same region of the phase diagram, namely for 
J/g = 4 (S = J/(2g) = 2 in our notation; tempera- 
tures in our paper are also rescaled by a factor 2 respect 
to those in Rastelli et al due to a different choice of the 
energy scales; see Ref4 for details). For that values of 
J / g Rastelli et all found a single peak in the specific heat 
and single minimum in the fourth order cumulant, whose 
finite size scaling is consistent with a first order phase 
transition, concluding that the system presents a unique 
first order phase transition between the (low tempera- 
ture) striped and the (high temperature) tetragonal liq- 
uid phases. On the other hand, we showed in Ref.— that 
for (5 = 2 the specific heat presents two distinct max- 
ima and the fourth order cumulant two distinct minima, 
consistent with the existence of an intermediate phase, 
characterized by orientational order and positional dis- 



order, evidenced in the finite size scaling behavior of the 
static structure factor. The existence of an intermediate, 
nematic phase, is consistent with one of the two possi- 
ble scenarios predicted by a continuum approximation 
for ultrathin magnetic films in Refs.— tt. Moreover, we 
showed that the finite size scaling of the different ther- 
modynamical quantities around the nematic-tetragonal 
phase transition are consistent with a first order tran- 
sition, while those observed around the stripe-nematic 
transition present unusual features, some of them re- 
sembling a Kosterlitz-Thouless one (like a saturation in 
the specific heat maximum and a continuous change in 
the internal energy) while others are also characteristic 
of a first order phase transition. In particular, one of 
those characteristics is the existence of diverging free en- 
ergy barriers between the nematic and the stripe phases, 
which generate a strong meta-stabilities when the sys- 
tem is cooled or heated through the transition tempera- 
ture. To correctly characterize the finite size scaling at 
this phase transition thermodynamical quantities must 
be averaged for every temperature over a large single MC 
run, thus allowing the system to sample properly the dif- 
ferent phases. Otherwise, as explained by Challa et al&, 
the first-order character of the transition introduces pro- 
nounced meta-stabilities with the system spending most 
of the time in one of the two phases (see Ref £ for details) . 
Hence, for every temperature we first let the system run 
over a transient period of Mi Monte Carlo Steps (MCS) 
and then we averaged the different quantities over a sin- 
gle MC run of M 2 MCS, with M x ranging between 10 6 
and 2 x 10 7 , and M 2 ranging between 2 x 10 7 and 2 x 10 8 ; 
we shall call this procedure protocol A. We believe that 
the absence of the nematic-stripe phase transition (and 
therefore the evidence of the nematic phase) in Rastelli et 
al paper is due to the particular MC protocol they used, 
which introduced a strong bias in their results. They took 
averages over eight independent MC runs of 10 5 MCS 
for each temperature, taking the initial configuration for 
every temperature as the final configuration of the pre- 
vious one, disregarding 10 4 MCS for thermalization 1 -; we 
shall call this procedure protocol B. This protocol is sim- 
ilar to a finite rate heating (cooling) procedure from low 
(high) temperature and the system can get trapped in 
the meta-stable striped (nematic) state, thus hiding the 
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FIG. 1: (Color on-line) Specific heat C, fourth order cumulant 
V and internal energy per spin u as a function of temperature 
for J/g = 4 (S = 2) and L = 48, obtained from different MC 
simulation protocols (see text for details). 



low temperature phase transition, whose associated free 
energy barriers are very high for system sizes used in that 
work. Indeed we checked this assumption by computing 
the moments of the energy, namely the internal energy 
per spin u, the specific heat C and the fourth order cumu- 
lant V as a function of the temperature (see definitions 
in Ref.-) using Rastelli et al MC protocol for 8 = 2 and 
L = 48, both heating from the low temperature stable 
phase and cooling from the high temperature one. The 
results were averaged over 8 (as in Rastelli et al paper) 
and 50 independent runs; no qualitative differences were 
observed. The obtained results for 50 independent runs 
are compared in FigQ]with the results obtained in Ref 4 



for the same parameters values. 

We see that the low temperature maximum of C and 
the minimum of V are almost undetectable by Rastelli 
et al MC protocol, the only noticeable effect being some 
enlarged fluctuations in the cooling procedure and a very 
small shoulder for the heating result at the left of the high 
temperature specific heat maximum (the inset of FigQ] 
remarks this effect) or cumulant minimum, but clearly 
separated from the actual low temperature transition. 
Moreover, the behavior of u shows clearly that the sys- 
tem gets trapped in a metastable super-cooled (super- 
heated) state, as in a finite rate cooling (heating) pro- 
cess (compare with Fig. 19 in Ref. 3 -). Indeed we observed 
that for some individual runs in the heating procedure 
the specific heat presents a small secondary maximum of 
the specific heat at low temperatures, but those maxima 
are a result of a spinodal instability of the super-heated 
state and therefore are strongly shifted to the right of the 
actual maximum. Moreover, those spinodal maxima ap- 
pear located at different random temperatures and they 
are smeared out when averaged, giving rise to the ob- 
served shoulder. 

Even when the analysis in Ref.— was carried out for 
the particular value 6 = 2, theoretical results^ lead us 
to expect the presence of the nematic phase to be found 
for a wide range of values of 8, although possibly for even 
narrower ranges of temperatures than for 8 = 2. In that 
case, the above described spurious effects should become 
worse. Indeed, the nematic phase and the associated be- 
havior in the specific heat and cumulant has been found 7 
using protocol A for other values of 8, including 8 = 2.25 
(J/g = 4.5), which corresponds to the h = 3 ground 
state region, as the J/g = 5 case analyzed by Rastelli et 
al (and close to it). 

Another related point concerns the scaling of the prin- 
cipal maximum of C for J/g — 5 in Rastelli et al work. 
They found that the maximum scales as L 12 instead 
of the L 2 scaling expected in a first order transition. A 
careful analysis of the orientational order parameter his- 
togram for 8 = 2 shows the presence of multiple peaks 
around this transition (See Fig. 7 in Ref.—) indicating the 
presence of a complex multiple phase structure in a very 
narrow range of temperatures; such structure is only de- 
tectable for large system sizes and can introduce a strong 
bias in the scaling properties for small sizes. The anoma- 
lous scaling observed by Rastelli et all is probably due to 
this finite size effect. 

The above considerations also apply to the case J/g = 
6. In this case the results of Rastelli et al (see Fig. 12 
in Ref fi) show rather clearly the presence of a secondary 
low temperature peak in the specific heat, whose height 
is independent of the system size, very similar to that 
observed for J/g = 4 (see Fig. 9 in Ref.—). There is also 
some indication of a vanishing secondary minimum in V, 
similar to that observed for J/g = 4. Although those 
anomalies are very small, they are probably affected also 
by the simulation protocol and by finite size effects, which 
in this system become stronger as J/g increases. Based 
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on the scaling of the specific heat maximum and on the 
apparently vanishing minimum of the energy cumulant 
the authors conclude that the transition is continuous 
for J/g > 6. However, such results can also be consis- 
tent with a weak first order transition that continuously 
fades out as J/g increases. If that would be the case, de- 
tecting the order of the transition becomes more delicate 
due to the above described finite size effects. Indeed, a 
Hartree approximation applied to the continuous version 
of Hamiltonian (HJ predicts a first order transition for a 
wide range of values of 8 J/g. Rastelli et al criticize this 
result claiming that such approximation wrongly predicts 
a first order transition also for g — 0, where it is known 
that it is second order. As is well known since the work 
by BrazovskiiS., that calculation breaks down for g = 0. 
Nevertheless, for a wide range of values of J/g > 1 the 
approximation gives excellent qualitative results. It is at 
least very suggestive that our numerical results seem to 
be in agreement with the predictions of this theory. To 
be more specific, at the transition a modulated solution 
of the form <p(r) = m cos (fcg • r) becomes stable and the 
amplitude of the equilibrium solution m(T) jumps dis- 
continously to a non zero value. Following Brazovskii^ it 
is possible to show that the amplitude at the transition 
is given by 

m c = m(T c ) oc y/r m (T c )/u, (2) 

where u is the coupling constant of the Landau energy 
and r m is the renormalized mass of the modulated solu- 
tion^. At the transition point r m oc (uk 2 ) 2 ^ 3 and 

,2/3 

m c oc -f- (3) 



In our model fco oc g/J and consequently the ampli- 
tude at the transition decreases continuously as g/J gets 
smaller within the range of validity of the approxima- 
tion. This shows that the transition becomes weaker as 
g/J diminishes. This is completely consistent with the 
difficulty in the numerical simulations to detect the first 
order nature of the transition for small values of g/J. 

Finally, for J/g = 1.7, Rastelli et al found evidence of 
a continuous phase transition with unusual critical expo- 
nents, from which they conclude that the h = 1 spin con- 
figurations (h being the width of the ground state stripes) 
undergo a continuous phase transition with a different 
universality class from the case J = 0. In Ref.— we pre- 
sented evidence that the system undergo a unique weak 
first order transition for J/g = 2 (S = 1); this value of the 
ratio also corresponds to the h = 1 region of the phase 
diagram (see phase diagram ini&) . These results are con- 
sistent with the presence of a second order transition line 
for small values of J / g that joins with continuous slope 
a first order transition line for larger values of J/g at a 
tricritical point in that region of the phase diagram (see 
phase diagram in 1 -). Hence, the unusual critical expo- 
nents observed by Rastelli et al are most probably due to 
a crossover effect related to the presence of that tricriti- 
cal point in the neighborhood of J/g = 1.7 and not to a 
new universality class transition line. 
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